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Abstract —Recently, hybrid-automatic-repeat-request (HARQ) 
systems have been favored in particular state-of-the-art commu¬ 
nications systems since they provide the practicality of error 
detections and corrections aligned with repeat-requests when 
needed at receivers. The qneueing characteristics of these systems 
have taken considerable focus since the current technology 
demands data transmissions with a minimum delay provision¬ 
ing. In this paper, we investigate the effects of physical layer 
characteristics on data link layer performance in a general 
class of HARQ systems. Constrncting a state transition model 
that combines quene activity at a transmitter and decoding 
efficiency at a receiver, we identify the probability of clearing 
the queue at the transmitter and the packet-loss probability at 
the receiver. We determine the effective capacity that yields the 
maximum feasible data arrival rate at the queue under quality-of- 
servlce constraints. In addition, we put forward non-asymptotic 
backlog and delay bounds. Finally, regarding three different 
HARQ protocols, namely Type-I HARQ, HARQ-chase combining 
(HARQ-CC) and HARQ-lncremental redundancy (HARQ-IR), 
we show the superiority of HARQ-IR in delay robustness over 
the others. However, we further observe that the performance gap 
between HARQ-CC and HARQ-IR is quite negligible in certain 
cases. The novelty of onr paper is a general cross-layer analysis 
of these systems, considering encoding/decoding in the physical 
layer and delay aspects in the data-link layer. 

I. Introduction 

Due to the dynamic nature of wireless media, transmission 
errors easily occur, especially when channel fading is too 
strong. Hence, in order to ensure an error-free (or error- 
minimized) transmission, error detection and correction tech¬ 
niques such as automatic-repeat-request (ARQ) and forward- 
error-correction (EEC) are engaged in certain classes of com¬ 
munications systems. Specifically in ARQ-furnished systems, 
the receiver feeds the transmitter with either a positive ac¬ 
knowledgment (ACK) or a negative ACK (NACK) after each 
data packet is received correctly or incorrectly so that the 
transmitter can either send the next packet or resend the same 
packet, respectively. On the other hand, in EEC each packet is 
encoded by introducing redundant data bits for the receiver 
to detect and correct errors. While EEC provides a fixed 
transmission rate, it does not adjust error correction capacity 
with respect to channel variations. However, ARQ can adapt to 
the channel variations by resending the transmitted packets. In 
that respect in ARQ systems, the transmission of a packet can 
be completed in a shorter period if the channel attenuation 
is weak, whereas it can be completed in a longer period if 
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the channel attenuation is strong. Thus, in order to provide 
better error correction performance, EEC and ARQ have been 
merged to form hybrid-ARQ (HARQ). 

Several HARQ protocols have been suggested and ana¬ 
lyzed. Eor instance, originally proposed by Wozencraft and 
Horstein HI, Type-I HARQ (HARQ-Tl) relies on a blended 
process of error detections and error corrections while the 
retransmission of one packet is initiated if this process fails. 
Subsequently, unlike HARQ-Tl in which the past erroneously 
acquired signals are discarded, Type-II HARQ was introduced 
in 0. Additional information is transmitted to the receiver 
when the receiver fails decoding the acquired signal. On that 
account, one of the leading studies introduced an adaptive- 
feedback coding strategy that performs Type-II HARQ with 
incremental redundancy (HARQ-IR) 13]. Likewise, a widely- 
known method of packet combining technique, named Type- 
II HARQ chase combining (HARQ-CC), has been analyzed 
in HI. In HARQ-CC, the repeated transmissions of a packet, 
before being decoded, are linearly added in order to obtain 
a power gain. Erom then on, Type-I and Type-II HARQ pro¬ 
tocols have been analyzed from various perspectives such as 
coding performance 0-Q, outage probability fSl, and power 
optimization a. Ho). Eurthermore, a link-level performance 
comparison of these two techniques has been investigated in 
M, and it was shown that HARQ-IR can be significantly 
better than HARQ-CC, especially for high channel-coding 
rates and high modulation orders. 

Meanwhile, since HARQ systems may require the retrans¬ 
missions of packets, delay-sensitive traffic concerns have be¬ 
come a research focus as well. Eor instance, considering adap¬ 
tive modulation techniques. Villa et al. investigated latency- 
constrained networks CD. In another line of research, the 
author in m formulated the design of an HARQ system as 
a stochastic control problem where the transmitter controls 
the code rate of each packet in order to achieve a minimum 
overall average delay. Similarly in m, effective capacity na, 
a dual of effective bandwidth m, has been characterized 
as the performance metric. A two-dimensional continuous¬ 
time Markov channel model has been employed in order to 
partition the instantaneous data rate received at the destination 
into a finite number of states each of which is representing a 
mode of operation of the HARQ scheme. In addition, Huang 
et al. conducted the analysis of queueing performance of 
HARQ systems by modeling the wireless channel as a two- 
state continuous Markov process and derived the queueing 
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Fig. 1. Channel model. The transmitter initially stores the data packets in its 
buffer, subsequently encodes, modulates and forwards each data packet to the 
receiver through the wireless channel. Then, the receiver feeds the transmitter 
with ACK/NACK depending on its decoding performance. 


delay from generating functions of the finite-buffer-capacity 
system ini. Using the framework of the statistical network 
calculus, the authors investigated the delay performance of 
block codes with more or less redundancy ifTSll . Finally, we 
note that taking the first order expansion of the effective 
capacity under loose quality-of-service (QoS) constraints into 
account, the authors in HU studied the impact of deadline 
constraints, outage probability, and QoS constraints on the 
energy efficiency in HARQ systems, and compared HARQ- 
CC and HARQ-IR protocols. 

Different than the above studies, we evaluate the queueing 
performance of a broad range of HARQ wireless systems 
by modeling the transmitter queue activity aligned with the 
decoding performance at the receiver as a discrete M-state 
Markov process with a transmission deadline strategy for 
each data packet. To the best of our knowledge, a cross-layer 
examination of a general class of HARQ systems, regarding 
the effect of wireless channel attenuation on decoding perfor¬ 
mance in physical layer together with queueing performance 
in data-link layer, has not been performed. Singly, we can 
express our contributions with the following; 

• We identify an approach and construct a state-transition 
model for any HARQ system to investigate how advanced 
transmission schemes at the physical layer affect the 
queueing performance at the data link layer. 

• Using the above model, we analyze the interplay between 
the packet-loss probability at the receiver due to transmis¬ 
sion deadlines, and the steady-state probability of clearing 
the transmitter queue. 

• We obtain a closed-form expression for the effective 
capacity of HARQ systems, which provides the maximum 
feasible arrival rate at the transmitter queue under QoS 
constraints such as the asymptotic probabilities of buffer 
overflow and delay violation. 

• Considering the stability condition of the queue, we 
analyze the non-asymptotic backlog and delay bounds of 
these systems. 

In this paper, we provide a general framework to investigate 
the interactions between the QoS concerns, and the coding 
and modulation techniques embedded in HARQ protocols. 


For instance, the adoption of linear coding and modulation 
schemes in HARQ protocols and their delay performance can 
be easily analyzed using this conduct. The organization of 
the rest of the paper is as follows. In Section [III we describe 
the system model. Then, we discuss in Section |III] the sys¬ 
tem performance regarding packet-loss probability, effective 
capacity, and non-asymptotic performance bounds. Finally, we 
provide the conclusion in Section |IV] and the relegated proofs 
in Appendices and |B] 

H. System Model 

In this section, we present in detail the transmission and 
channel model, a general view of HARQ systems, and the 
state transition model of these systems. 

A. Transmission and Channel Model 

We consider a point-to-point transmission system in which 
one transmitter and one receiver communicate over a wireless 
fading channel as seen in Figure [T| In this model, it is assumed 
that the transmitter initially divides the available data into 
packets of n bits and stores them in its buffer. Then, it performs 
the encoding, modulation, and transmission of each packet in 
time-slots of T seconds with the first-come first-served policy. 
It is further assumed that each packet should be transmitted 
in maximum MT seconds due to a transmission deadline for 
each packet where M is an integer. Specifically, if a packet is 
not received correctly by the receiver at the end of the 
slot, it will be discarded from the transmitter queue. During the 
transfer of each packet, if the packet is successfully decoded at 
the end of any slot, the receiver feeds the transmitter with an 
ACK, and the transmitter removes the packet from its queue. 
Otherwise, the receiver sends a NACK, and the transmitter, 
depending on the transmission protocol, either repeats the 
same encoded packet or sends the next partition of a larger 
codeword in the next slot. We note that a removal of any 
data packet from the buffer, either as a result of a successful 
decoding by the receiver or due to the transmission deadline, 
is considered to be a packet service from the queue. 

All along, the discrete-time channel input-output relation in 
the symbol instant is given as 

y{t) = h{t)x{t) + w{t) for f = 0,1, • • • , 

where x{t) is the complex channel input and y{t) is the 
complex channel output. Above, {^(t)} forms an independent 
and identically distributed (i.i.d.) sequence of additive zero- 
mean circularly symmetric, complex Gaussian random noise 
variables with a variance E{|w(f)p} = cr^. Furthermore, 
h{t) denotes the fading coefficient between the transmitter 
and the receiver with an arbitrary marginal distribution and 
a finite variance E{|/i(f)p} = E{ 2 (f)} = < oo. Note 

that, here and throughout the paper, z{t) = denotes 

the magnitude-square of the fading coefficients. We further 
assume that channel side information (CSI) is available at the 
receiver, i.e., the receiver knows the actual value of h{t), while 
the transmitter is only aware of channel statistics. Besides, due 
to a limited power budget, the channel input is subject to the 
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Fig. 2. State-transition model with M states. The number of bits removed from the transmitter queue is n in State 0, while it is 0 in other states. 


following average power constraint: E{|a;(f)p} < P/B where 
B is the available channel bandwidth. Since we assume that 
B complex symbols per second are transmitted, the average 
power of the system is constrained by P. Finally, we consider 
a block-fading channel model and assume that the fading 
coefficient stays constant for a slot duration of T seconds 
and changes independently from one slot to another, i.e., 
h{lTB) = h{lTB + 1) = ... = h{{l -t- l)TB -1) = hi 
in the slot, and = z;. This channel model has been 
widely used to represent slowly-varying, flat fading channels 
II 20 I and references therein]. 

In the above model, the instantaneous ergodic channel 
capacity in the slot is given by 

Cl = i?log 2 {1 + 7 Zi} bits/sec, for / = 0,1, • • • , (1) 

where the average transmitted signal-to-noise ratio (SNR) is 
7 = -g^. In ([T]), the capacity is achievable when x{t) is zero- 
mean Gaussian-distributecQ, i.e., x{t) CAf{0, P/B) El\ Ch. 
9.1]. Since a block-fading channel is considered, we denote the 
channel fading power and the instantaneous channel capacity 
in the slot by zi and Ci, respectively. Additionally, we 
remark that if the transmitter would know the instantaneous 
channel fading coefficient hi, it could send the data with a rate 
equal to the instantaneous channel capacity. Since in practice 
hi is generally not available at the transmitter, we consider a 
more realistic scenario, where a wireless link is equipped with 
an HARQ protocol, which can be modeled as a discrete-time 
finite state Markov chain. 

B. HARQ 

In general, we assume that before being delivered into 
the channel, a packet of n bits is encoded and modulated 
into a codeword of length either TB or MTB complex 
symbols. We can easily deduce that since one fading slot 

'Codewords can be composed of symbols from any modulation technique 
with the input distribution Q{x). Then, the capacity in (T) is replaced with 
the achievable rate expressed as the mutual information between x and y, 
i.e., I{x\y). The physical layer framework of this study and the analysis in 
the rest of the paper can be easily adopted to any practical linear modulation 
technique by using Q{x) and I{x\y). 


is T seconds, and the available bandwidth is B Hz, a total 
of TB complex symbols can be transmitted in one slot. 
Therefore, a codeword of TB complex symbols is assumed to 
be repeatedly transmitted in each slot until either it is decoded 
by the receiver or the transmission deadline is reached. HARQ- 
T1 and HARQ-CC are two of these protocols that perform 
iterative transmissions of data packets. On the other hand, a 
codeword of MTB complex symbols is abstracted into M 
sub-codewords that are transmitted consecutively in each slot 
again until either it is decoded or the transmission deadline is 
reached, e.g., HARQ-IR. 

Particularly, the first partition (codeword or sub-codeword 
of length TB complex symbols) is sent to the receiver in 
the first attempt If a successful decoding is performed, the 
receiver sends an ACK to the transmitter, and then the packet 
is removed from the queue. On the other hand, if a decoding 
failure occurs, the receiver sends a NACK to the transmitter 
and requests the second partition. Then, the transmitter sends 
the second partition in the second attempt. Equivalently, fol¬ 
lowing a successful decoding in any attempt, the packet is 
removed from the queue; otherwise with a decoding failure, 
the corresponding next partition is transmitted in the next 
attempt. However, at the end of the attempt, the packet 
is removed from the queue without waiting for any ACK or 
NACK since the transmission deadline for a packet is reached. 
We finally note that due to a decoding failure at the end of the 
attempt, if a packet is removed from the queue and not 
transmitted again, we deem it as a lost packet at the receiver 
in order to differentiate it from the other packet removals. So, 
the ratio of the lost packets to the total number of packets 
served by the transmitter is a very critical QoS measure since 
it entails the average reliable transmission rate to the receiver. 

C. State Transition Model 

For an analytical presentation, we model the queue activity 
at the end of each slot as a discrete-time Markov process. As 
seen in Fig. |2] following a decoding success at the end of 

^In order to differentiate any slot from the slot in which a packet is 
transmitted for the first time, we use attempt rather than slot, unless otherwise 
needed for clarity. 
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the attempt, the system enters State 0 with probability 1- 
Pm-i, whereas if a decoding failure occurs the system enters 
State m with probability Pm-i where m G {I,-- - , M — 1}. 
On the other hand, regardless of the decoding result at the end 
of attempt, the system goes to State 0 with probability 1 
since the transmitted packet is removed from the queue due to 
the transmission deadline. We note that a packet transmitted in 
the attempt is decoded either correctly with probability 
1-pM-i or incorrectly with probability pm-i- Accordingly, 
we have the following state transition matrix P: 


end of the slot for 1 < m < M is expressed as Co) 


B 




i=0 


bits/sec. 


Above, while Cj is the channel capacity in the first slot, 
can be considered as the average channel capacity in 
the first k slots. Noting that the receiver can decode a packet 
at the end of slot when < C;we have the 
following state transition probabilities: 
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Now, noting that the analysis in the paper is for a general 
HARQ protocol, we consider the following three different 
HARQ protocols in a RayleiglU channel fading environment 
for a profound examination. 

I) Type-I HARQ (HARQ-Tl): In this protocol, the transmit¬ 
ter encodes and modulates its data packets, of n bits, into one 
codeword by using a channel code of a codebook C G 
of length TB over the complex numbers. The codeword is 
transmitted over the channel in one slot. We assume that if the 
transmission rate is lower than or equal to the instantaneous 
channel capacity in O, i.e., ^ < Cl, the receiver is able 
to decode the codeword reliably. On the other hand, if the 
transmission rate is greater than the capacity, i.e., ^ > Ci, 
the receiver cannot recover the packet correctly. Featuring a 
receiver that does not perform accumulation of data, i.e., that 
uses only the current received signal in any slot in order to 
retrieve the packet, the decoding failure probability in any slot 
is given by 

p = Pr{^>G} = Pr{K >zi} = l- e-C-l 

where the probability density function of exponentially- 
distributed zi (i.e., Rayleigh distributed hi) is fzi{zi) = 
and K = (2to — l)/ 7 . Since any state transition does 
not depend on the past transitions, we have p = po = ■ ■ ■ = 




Pr < 

1 n 

[(m-l)T ' 

Pr j 



\-m — 2 


} 


7 I m 




where m G {2, ■ ■ ■ ,M}, and po = 1 — Above, 7 ( 0 , b) 

is the lower incomplete gamma function where 7 ( 0 , b) = 
Jq and k is as defined in Section riTCll Note 

that the sum of exponentially-distributed random variables is 
gamma-distributed. 

3) Type-II HARQ-IR: In this protocol, different than the 
other two protocols, we assume that the transmitter encodes 
and modulates its data packets, of n bits each, by using a 
channel code of a codebook C G of length MTB over 

the complex numbers. Then, the transmitter divides the code¬ 
words into M partitions of the same length with TB complex 
symbols. In each slot, one partition is sent to the receiver, and 
the receiver decodes the message using the current partition 
combined with the previously received partitions related to the 
encoded data packet. Hence, following the result in 12^ . we 
can see that the receiver can decode the transmitted data packet 
at the end of the slot, if y < ^ (1 + 

Subsequently, we can express the state transition probabilities 
for m G { 2 , • • • , M} as follows: 


Pr| 


Pm — 1 — 
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■ 1Z1+ 
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Pr I 7 K -b 1 > (1 + 


Pm-i- 

2) Type-II HARQ-CC: In this protocol, we again assume 
that the transmitter encodes and modulates its data packets, of 
n bits each, by using a channel code of a codebook C G 
of length TB over the complex numbers. However, different 
than the receiver structure in Section lEI] we devise a 
receiver that makes use of the received signals in the earlier 
slots by implementing maximum ratio combining. Hence, the 
average accumulated mutual information at the receiver at the 

^The Rayleigh fading channel as sum ption fits well with HARQ schemes in 
practical mobile wireless systems Qo). However, other channel models can 
be adopted easily into this study. 


where po = 1 — e and k is as defined in Section Hl-C II 

HI. System Performance 

In this section, we concentrate on the performance measures 
of the HARQ systems with the above defined state-transition 
model. We initially obtain the packet-loss probability at the 
receiver, and then determine the effective capacity which 
provides the maximum sustainable arrival rate at the buffer 
under certain QoS constraints. We finally provide the non- 
asymptotic performance measures in the form of backlog and 
delay bounds. 
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Fig. 3. Reliable throughput, bits/sec., v.s. number of bits 

encoded and modulated into TB symbols in one slot, n, when 7 = 5 dB. 


Fig. 4. Steady state probability of State 0, tto, as a function of transmissions 
deadline, M, for three different protocols when 7 = 0 and 10 dB. 


A. Packet-Loss Probability 

Let TT = [tto,--- jTTm-i]^ be the steady-state probability 
vector of the above state transition model in Fig. |2] where 
TTj = 1 and TT = Ptt. Thus, we can obtain the steady- 
state probabilities as follows: 


i-l 

TTj = ttq Pj and ttq 
i=o 


-1 I t^M—1 t—tz—1 

1 “b Y[j=oPj 


where i G {I,-- - ,M — 1}. Recalling that at the end of the 
slot, data packets are removed from the queue due to a 
decoding failure with probability pm-i and due to a decoding 
success with probability 1 —pM-i, we can easily express the 
packet-loss probability, i.e., the ratio of the lost packets to the 
total number of packets that are served by the transmitter, with 


Plost — 


TTo 


M-2 M-1 

PM-i 'Wpj= n pp 


Note that while ttq is the steady-state probability of removing 
a packet from the queue in one slot, piost is the probability of 
a packet not reaching the receiver. 

We employ the aforementioned HARQ protocols, and plot 
the average reliable throughput that arrives at the receiver, i.e., 
bits/sec, as a function of n, i.e., the number of 
bits that are encoded and modulated into TB symbols in one 
slot, when 7 = 5 dB in Fig. |3] We consider the following 
settings: The average channel fading power is = 1, the slot 
duration is T = 100 psec, and the transmission bandwidth is 
B = I MHz. We observe that while HARQ-IR significantly 
outperforms the other two protocols, HARQ-Tl is the least 
achieving protocol for all n. Fixing a constant number of 
bits served in all protocolfl we further compare the queue 
clearing performances of all the protocols, ttq, as a function 
of the transmission deadline, M, when 7 = 0 and 10 dB in 
Fig. a We can see for any 7 that the probability ttq decreases 
generally with increasing M, then it converges after certain M 


'^Unless stated otherwise, here and throughout the paper, we consider 
the given transmission settings, and choose n that maximizes the reliable 
throughput at the receiver when FIARQ-Tl is employed. 



Fig. 5. Packet-loss probability, piost, as a function of transmissions deadline, 
M, for three different protocols when 7 = 0 and 10 dB. 


in all HARQ protocols. All along, HARQ-IR converges very 
quickly, while HARQ-CC performs very close to HARQ-IR. 
On the other hand, HARQ-TI accomplishes considerably less 
in general compared to the Type-II HARQ protocols for all 7 
values. Finally, we plot the packet-loss probability as a func¬ 
tion of M in Fig.|5] Unlike ttq, piost is continuously decreasing 
with increasing M. HARQ-IR has a significantly lower packet- 
loss probability, and the performance gap between HARQ-IR 
and the others increases dramatically with again increasing M. 


B. Effective Capacity 

Even though it seems to be advantageous to increase M in 
order to decrease piost, we can promptly infer that the average 
delay of a packet in the buffer will increase with increasing M. 
Therefore, we propose the effective capacity that characterizes 
the asymptotic decay rate of buffer occupancy. It identifies the 
maximum constant arrival rate that a given service process can 
support in order to guarantee a desired statistical QoS specified 
with the QoS exponent 0 1151. Defining Q{t) as the stationary 
queue length at time t, and 9 as the decay rate of the tail 
distribution of the queue length Q{t), we can express the fol¬ 
lowing: lim^^oo *°s _ _() Therefore, we have the 

following approximation for larger q: Pi{Q{t) > q) ^ 
which means that larger 9 refers to strict buffer constraints. 
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Fig. 6 . Effective capacity, ps{S), as a function of decay rate of the tail 
distribution, B, for three different protocols when 7 = 0 dB. 


Fig. 7. Effective capacity, ps{S), as a function of decay rate of the tail 
distribution, 6, for three different protocols when 7 = 5 dB. 


and smaller 9 implies looser constraints. Furthermore, it is 
shown in ll^ that Pr(D(t) > d) < c^/Fv{Q{t) > q) for 
constant arrival rates, where D{t) denotes the steady-state 
delay experienced in the buffer, and c is a positive constant. 
In the above formulation, we have q = ad. Therefore, 
effective capacity provides us with the maximum arrival rate 
when the system is subject to the statistical queue length or 
delay constraints in the forms of Pr((5(f) > q) < or 

Pi{D{t) > d) < respectively. For a given QoS 

exponent 6 > 0, effective capacity, ps{d), is given by 

= = ( 2 ) 

where S(t, t) = J2\=t+i time-accumulated service 

process, and r{k) for k = 1,2, ... is the discrete-time, station¬ 
ary and ergodic stochastic service process. We note that A{6) is 
the asymptotic log-moment generating function of S'( 0 , t), and 
is given by A{9) = limi_i.oo 7 ^ogE Henceforth in 

the next result, we provide the effective capacity for any given 
HARQ protocol. 

Theorem 1: For the HARQ system with the state-transition 
model given in Section ITl-Cl the effective capacity for a given 
QoS exponent 9 is given by 

Ps(fi') =-^log (poe“”®y*) bits/sec (3) 


where y* is the only unique real positive root of f{y) where 

M PM-2---P1 


f{y) = 2 /“ - 

Po 


pM-le-iM-l)ne 


2=1 


(4) 


— XTtO 

Po^ 


-y 


Proof: See Appendix lAl 

The real positive root of f{y) can be found by using 
numerical techniques. For instance, bisection method can be 
efficiently used to find the solution since it has only one real 
positive root. 

Employing, again, the above HARQ protocols, we plot the 
effective capacity, ps{d), as a function of the decay rate, 
9, when 7 = 0 dB and 7 = 5 dB in Fig. | 6 ] and Fig. 


|7] respectively. Above, we set the transmission deadline to 
M = 4. In both figures, the superiority of HARQ-IR is clearly 
seen, and HARQ-CC performs close to HARQ-IR. We further 
see that while the effective capacity of each protocol goes 
to the average service rate of the same protocol with 
decreasing 9, all of them go to the same minimum service 
rate with increasing 9. 


C. Non-asymptotic Performance Bounds 

So far, we have considered the effective capacity that deals 
with asymptotic characteristics. However, non-asymptotic per¬ 
formance bounds regarding the statistical characterizations of 
backlog and delay are also of interest for system designers. 
For a non-asymptotic analysis we use the framework of the 
stochastic network calculus lED-iIlTl. Following the model in 
II 24 I Definition 7.2.1], we define a statistical affine bound for 
the above channel model for any 9 as follows; 

g-6»[ps(6()(t-T)-(Ts(6()] 

where ps{9) is the effective capacity and crs(9) is a slack term 
that defines an initial service delay. Although the expression 
in 0 seems to be an upper bound, due to —9, in fact 0 
is a lower bound on the expected amount of service. Now, 
using Chernoff’s lower bound Pr{X < x} < 
for 9 > 0 , the so-called exponentially bounded fluctuation 
model described in ll28l with parameters ps{9) > 0, 6 > 0, 
and Pr {S'(r,f) < ps{9){t — t) — b} < e{b) follows, where 
e{b) = g®‘^s(e)g-eb specifies an exponentially decaying deficit 
profile of the service. A sample path guarantee follows by the 
use of the union bound as: 


Pr{3r e [0,f] : < Ps(6»)(f - t) - b} < efb) 


where 


e'{b) 


1 — 


e 


-9b 


( 6 ) 


and p'g{9) = ps{9) — 5 with free parameter 0 < ^ < ps{9)-a 
for an arrival rate a. For a detailed derivation, we refer the 
interested reader to Wl- All in all, perceive that the backlog 
Q{f) = maxT-g[Q q{a(f — r) — S{t, t)} has a statistical bound 
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Fig. 8. Delay bound, d, as a function of the violation probability, s', for 
three different protocols when 7 = 0 and a = 0.41 Mbits/sec with n = 82. 



Fig. 9. Delay bound, d, as a function of the violation probability, e', for 
three different protocols when 7 = 5 and a = 0.81 Mbits/sec and n = 155. 


Fig. 10. Delay bound, d, as a function of signal-to-noise ratio, 7, when 
M = 4 and a = 0.16 Mbits/sec. 



Fig. 11. Delay bound, d, as a function of signal-to-noise ratio, 7, when 
M = 3 and a = 0.18 Mbits/sec. 


q = max^g[o.t]{a(i - t) - [p's{d){,t - r) - b]+} that may 
fail with probability Pr{Q(f) > g} < s'{b), where [a;]+ = 0 
if a; < 0, and [a;]+ = x otherwise accounts for the fact that 
5'(t, f) > 0, generally. If a < p's{d) for stability, 


follows for all t. Accordingly, the delay bound Pr{D(f) > 
d} < £'{b) can be expressed with d = q/a. In (|7]i, b/{ps{d) — 
(5) provides us with the initial latency caused by the variability 
of the service. Finally, by inversion of (l6]l, b can be easily 
obtained for any given e' by 

b = crs{0) - i (log(£') + log (1 - e~^^)) . (8) 

As for the existence of the slack term as {9) in ([8]l, consid¬ 
ering the result ll24l 7.2.6.ii], we provide the following lemma: 

Lemma 1: If S{T,t) has an envelope rate ps{d) < oo, for 
every e > 0, there exists as{9) < oo such that S{T,t) is 
{as{d), ps{9) — e)-upper constrained Q. 

Proof: See Appendix iBl 

As seen in Fig. [8] and Fig. |9] we plot the delay bound 


threshold, d, as a function of the violation probability, e', when 
7 = 0 dB and 7 = 5 dB when the arrival rates are a = 0.41 
Mbits/sec and a = 0.81 Mbits/sec while n = 82 bits/slot and 
n = 155 bits/slot are encoded and modulated into 100 sym¬ 
bols, respectively. We observe that HARQ-Tl has a very poor 
performance when compared with the other Type-II HARQ 
protocols. For instance, d is around 78 msec when HARQ- 
Tl is utilized, whereas it is around 5 msec when the Type-II 
HARQ protocols are engaged for e = 10“® and 7 = 0 dB. 
Similarly, when 7 = 5 dB, d is reduced from around 11 msec 
to around 4 msec. We remark that the arrival rates are different 
and are close to the average service rates in both plots. We 
again observe that HARQ-CC and HARQ-IR achieve delay 
bounds very close to each other. Furthermore, setting a fixed 
arrival rate a = 0.16 Mbits/sec and a = 0.18 Mbits/sec, we 
plot the delay bound threshold, d, as a function of SNR, 7, 
in Fig. [To] and Fig. [TT] when M = 4 and 3, respectively, 
while e' = 10“® and n = 36 bits/slot. As expected, the delay 
bound decreases with increasing 7. However, the decrease in 
d converges in all HARQ protocols. This is due to the fixed 
transmission slot length, T. It is worth mentioning that when 
the delay is the main concern, expending transmission energy 
above some value is useless, and that when there is enough 
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Fig. 12. Delay bound, d, as a function of data arrival rate, a, when 7 = 5 
dB. 


Fig. 13. Delay bound, d, as a function of data arrival rate, a, when 7 = 10 
dB. 


energy provided for data transmission, employing a simplified 
technology is more profitable for system designers. As for the 
delay performance regarding the variations in data arrival rates 
at the buffer, we plot d as a function of a in Fig. [12] and Fig. 
fm when 7 = 5 dB and 10 dB while n = 155 bits/slot and 
252 bits/slot, respectively. We again note e' = 10“®. In all 
protocols and cases, when the arrival rate is greater than the 
average service rate, the delay threshold goes to infinity since 
the system is not stable any more. The performance of HARQ- 
T1 is noticeably degraded with increasing a when compared 
to other HARQ protocols. However, a system furnished with 
HARQ-Tl is more appreciated when arrival rates are lower. 

IV. Conclusion 

In this paper, we have considered HARQ systems under 
QoS constraints such as asymptotic and non-asymptotic delay 
and backlog bounds. We have addressed a practical setting in 
which data packet transmission under a transmission deadline 
is performed. We have initially established a state-transition 
model to define the queueing characteristics of the system, 
then we have identified the steady-state probability of clearing 
the queue at the transmitter and the packet-loss probability at 
the receiver. We have further derived the effective capacity that 
features the arrival rate at the transmitter queue under QoS 
constraints. We have then furnished the analysis with the non- 
asymptotic delay and backlog bounds with respect to a given 
violation probability. We have finally presented the numerical 
results comparing three different HARQ protocols. 

Appendix 

A. Proof of Theorem Q] 

In ||24| Chap. 7, Example 7.2.7], it is shown for Markov 
modulated processes that 

^ = ^logsp{P<?!>(6»)} = ilogsp{T} 

where sp{T} is the spectral radius of the matrix T, P 
is the transition matrix of the underlying Markov process, 
and (j){9) = diag{7(^*)! is a diagonal matrix. 


components of which are the moment generating functions of 
the processes in M states. The rates (number of bits leaving 
the queue) supported by the above channel model with the 
state transition model described in the previous section can 
be seen as a Markov modulated process and hence the setup 
considered in ll24ll can be applied immediately into our setting. 
We have n bits served from the queue in State 0 while 
we have 0 bit served in other states. Therefore, we have 


(j){6) = diagje 


—n9 


, 1 ,. 


T = 


((1 -Fo)e 
Poe-"® 
0 


— n9 


, 1}. Then, we have 

1-pi l-pM-2 l\ 

0 00 

Pi ■■■ 0 0 


(9) 


V 0 0 ■ ■ ■ Pm-2 oJ 

Note that T in (|9]) is a non-negative matrix, i.e., T > 0, with 
each element greater or equal to zero, i.e., Vij > 0, where 
T = [vij]. f2^ Chap. 8, Corollary 8.3.3] states that if T > 0, 


= [xq, ■ ■ ■ ,Xm-i] > 0, and x 7^ 0 then 


sp(T)=max min — 

x>0 0<i<M-l Xi 

min|(l-po)e-"® 


^ M-1 


= max 

x>0 
x^O L 


XM-l + - Pi)^i 

Xo 


—nO Xq Xi X]if—2 

Poe — ,pi —, • • • , pm-2 - 

Xl X2 XM-1 


( 10 ) 


Now, let us assume that the maximum of the minimum in (ITOi) 

is obtained whenever 

—nd ^0 X^—2 

Poe — =pi — = ---=pm-2 - 

Xl X2 XM-1 


XM-1 + ^(1 “ Pi)^ 


= (l-po)e-"® + 


( 11 ) 


Xo 


Note that whenever Xi for i = 0, • • • , M — 1 changes, the 
minimum in (ITOl i will decrease, but it will not be the maximum 
of the minimum values. Therefore, sp(T) will be obtained 
































whenever (fTTl i is provided. Holding the equality in (fTTl i. we 
have 

Pi-lP^-2- ■ -pi f XiV ^ 

Xi = xo — for 1=2 -•-M-l. (12) 

pI \xoJ 

Then, defining y = —, and inserting y and (fT2l i into the 

I— 

equality in (fTTTi . we have sp(T) = pQe~™y and 


M-2 


y =- 


l-PO , - Pi)Pi-l ■ ■ ■ Pi 1 


Po 


E 


PM-2 • • - Pl 


yy — 1/7X0 

Po^ 

1 


(13) 


M-l^-(M-l)ne ■ 


Po 


Reordering (fOl l. we obtain the expression f{y) in (|4]i. In order 
to analyze the roots of fijj), we invoke the following theorem: 


Theorem 2 (Cauchy’s Theorem): Let g{x) = x^ — 
— ■ ■ ■ — bn, where all the numbers bi are non-negative 
and at least one of them is non-zero. The polynomial g{x) has 
a unique positive root x* and the absolute values of the other 
roots do not exceed x* ISO). 

Note that f{y) has coefficients that are non-negative and 
at least one of them is non-zero. Therefore, there is only one 
unique real positive root of f{y), denoted by y*, which gives 
the spectral radius of T. 


B. Proof of Lemma Q] 

We can note from dU that for every e > 0, there exists 
to < oo such that 

i sup log < ti-ps{0) + e) 

^ T>0 


<h{-ps{0)Pe) (14) 

when t > to. Meanwhile, for all t < to, we can write 


isuplogL;e-®'^(^’^+‘) < 0. 

P r>0 


Let as{0) = to{ps{0) — e) such that for t < to 
EuplogL;e-®^("’"+‘) <O<f(-ps(0) + e) + as(0). (15) 

tf r>0 

From di and (fTSl) . we can easily see that S(t,t -f t) is 
{as{d), ps{0) — e)-upper constrained, i.e., for all 0 < f 

suplogL;e-®''(">"+‘) < -e[t{ps{9) - e) - as{9)]. 

T>0 


References 

[1] J. M. Wozencraft and M. Horstein, “Coding for two-way channels,” 
DTIC Document, Tech. Rep., 1961. 

[2] S. Lin and P. S. Yu, “A hybrid ai'q scheme with parity retransmission 
for eiTor control of satellite channels,” IEEE Trans. Commun., vol. 30, 
no. 7, pp. 1701-1719, 1982. 

[3] D. M. Mandelbaum, “An adaptive-feedback coding scheme using incre¬ 
mental redundancy (corresp.),” IEEE Trans. Inf. Theory, vol. 20, no. 3, 
pp. 388-389, 1974. 

[4] D. Chase, “Code combining-a maximum-likelihood decoding approach 
for combining an arbitrary number of noisy packets,” IEEE Trans. 
Commun., vol. 33, no. 5, pp. 385-393, 1985. 


[5] J.-F. Cheng, “Coding performance of hybrid arq schemes,” IEEE Trans. 
Commun., vol. 54, no. 6, pp. 1017-1029, 2006. 

[6] A. Avudainayagam, J. M. Shea, and A. Roongta, “Improving the 
efficiency of reliability-based hybrid-arq with convolutional codes,” in 
IEEE Military Commun. Conf (MILCOM). IEEE, 2005, pp. 448^54. 

[7] E. Malkamaki and H. Leib, “Performance of truncated type-ii hybrid arq 
schemes with noisy feedback over block fading channels,” IEEE Trans. 
Commun., vol. 48, no. 9, pp. 1477-1487, 2000. 

[8] P. Wu and N. Jindal, “Performance of hybrid-arq in block-fading 
channels: A fixed outage probability analysis,” IEEE Trans. Commun., 
vol. 58, no. 4, pp. 1129-1141, 2010. 

[9] D. V. Djonin, A. K. Karmokar, and V. K. Bhargava, “Joint rate and power 
adaptation for type-i hybrid ai'q systems over correlated fading channels 
under different buffer-cost constraints,” IEEE Trans. Veh. Technol., 
vol. 57, no. 1, pp. 421^35, 2008. 

[10] T. V. Chaitanya and E. G. Larsson, “Optimal power allocation for hybrid 
arq with chase combining in iid rayleigh fading channels,” IEEE Trans. 
Commun., vol. 61, no. 5, pp. 1835-1846, 2013. 

[11] P. Frenger, S. Parkvall, and E. Dahlman, “Performance comparison of 
harq with chase combining and incremental redundancy for hsdpa,” in 
IEEE Veh. Technol. Conf Fall (VTC-FALL), vol. 3. IEEE, 2001, pp. 
1829-1833. 

[12] T. Villa, R. Merz, R. Knopp, and U. Takyar, “Adaptive modulation and 
coding with hybrid-arq for latency-constrained networks,” in European 
Wireless Conf VDE, 2012, pp. 1-8. 

[13] A. Anastasopoulos, “Delay-optimal hybrid arq protocol design for 
channels and receivers with memory as a stochastic control problem,” 
in IEEE Int. Commun. Conf. (ICC). IEEE, 2008, pp. 3536-3541. 

[14] N. Gunaseelan, L. Liu, J. Chambeiiand, and G. H. Huff, “Performance 
analysis of wireless hybrid-arq systems with delay-sensitive traffic,” 
IEEE Trans. Commun., vol. 58, no. 4, pp. 1262-1272, 2010. 

[15] D. Wu and R. Negi, “Effective capacity: a wireless link model for 
support of quality of service,” IEEE Trans. Wireless Commun., vol. 2, 
no. 4, pp. 630-643, 2003. 

[16] C.-S. Chang and J. A. Thomas, “Effective bandwidth in high-speed 
digital networks,” IEEE J. Sel. Areas Commun., vol. 13, no. 6, pp. 1091- 
1100, 1995. 

[17] L. Huang and T. T. Lee, “Queueing behavior of hybrid arq wireless 
system with finite buffer capacity,” in IEEE Wireless Optical Commun. 
Conf (WOCC). IEEE, 2012, pp. 32-36. 

[18] R. Lubben and M. Fidler, “On the delay performance of block codes 
for discrete memoryless channels with feedback,” in IEEE Sarnojf Symp. 
(SARNOFF). IEEE, 2012, pp. 1-6. 

[19] Y. Li, G. Ozcan, M. C. Gursoy, and S. Velipasalar, “Energy efficiency 
of hybrid-ai'q systems under qos constraints,” in IEEE Inform. Sci. Syst. 
(CISS). IEEE, 2014, pp. 1-6. 

[20] R. A. Berry and R. G. Gallager, “Communication over fading channels 
with delay constraints,” IEEE Trans. Inf. Theory, vol. 48, no. 5, pp. 
1135-1149, 2002. 

[21] T. M. Cover and J. A. Thomas, Elements of information theory. John 
Wiley & Sons, 2012. 

[22] G. Caire and D. Tuninetti, “The throughput of hybrid-arq protocols for 
the gaussian collision channel,” IEEE Trans. Inf. Theory, vol. 47, no. 5, 
pp. 1971-1988, 2001. 

[23] L. Liu and J.-F. Chamberland, “On the effective capacities of multiple- 
antenna gaussian channels,” in IEEE Int. Symp. Inform. Theory (ISIT). 
IEEE, 2008, pp. 2583-2587. 

[24] C.-S. Chang, Performance guarantees in communication networks. 
Springer Science & Business Media, 2000. 

[25] F. Ciucu, A. Burchai'd, and J. Liebeherr, “Scaling properties of statistical 
end-to-end bounds in the network calculus,” IEEE Trans. Inf. Theory, 
vol. 52, no. 6, pp. 2300-2312, 2006. 

[26] Y. Jiang and Y. Liu, Stochastic network calculus. Springer, 2008, vol. 1. 

[27] M. Fidler and A. Rizk, “A guide to the stochastic network calculus,” 
IEEE Commun. Sitn’eys Tuts., vol. PP, no. 99, pp. 1-1, 2014. 

[28] K. Lee, “Performance bounds in communication networks with variable- 
rate links,” in ACM SIGCOMM Comput. Commun. Review, vol. 25, 
no. 4. ACM, 1995, pp. 126-136. 

[29] R. A. Horn and C. R. Johnson, Matrix analysis. Cambridge university 
press, 2012. 

[30] V. V. Prasolov, Polynomials. Springer Science & Business Media, 2009, 
vol. 11. 


9 







